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Abstract 

The perturbation theory is developed for joint statistics of the advanced and retarded 
Green's functions of the ID Schrodinger equation with a piecewise- const ant random po- 
tential. Using this method, analytical expressions are obtained for spectral dependence of 
the degree of localization and for the limiting (at t — )• oo) probability to find the particle at 
the point it was located at t = (Andeson criterion). Definition of the localization length 
is introduced. The computer experiments confirming correctness of the calculations are 
described. 

1 Introduction, Formulation of the Problem, and Main 
Results 

Mathematical problems arising in the physics of solid-state random systems are characterized 
by a complexity and absence of universal methods of analysis. In searching for such methods, 
an important role is played by strongly simplified models of disordered systems, with the ID 
single-particle ones being the most important among them. The heuristic significance of the ID 
models is, however, not the only one. The physical systems like J-aggregates, quantum wells, 
optical fibers, Bragg layered structures, etc. can be directly described by the ID models. In 
such systems, one may expect strong effects of disorder, which makes studying of the disordered 
ID models especially topical. 

In the theory of the simplest solid-state disordered models, one can distinguish the con- 
tinuous and discrete models[T]. The continuous models employ the Schrodinger equation 



[-<P/dx^ + U{x)]^ = Eip with one or another potential U{x), while the discrete models use 
random matrices of a model Hamiltonian. In spite of similarity of these two models, each 
of them has its own specificity. For example, when analyzing vibrations of disordered chains 
[21 [3] , the discrete model is used, whereas when studying propagation of electromagnetic waves 
in disordered layered structures, the continuous model looks more convincing. 

For the ID models, in a number of cases, mathematically correct methods of theoretical 
analysis can be proposed [U E], [1]. Of particular interest are the cases, when one and the 
same method appears to be suitable for analysis of several different model problems, and the 
method proves to be, to a certain extent, universal. It is noteworthy, in this connection, that 
the perturbation theory for the joint statistics of the advanced and retarded Green's functions 
used in this paper for analysis of the continuous disordered model, has been successfully used 
previously for studying the discrete models [6l [7j . 

Let us pass to the problem studied in this paper. Consider a ID continuous disordered 
model with a peacewise-constant random potential U{x) equal to m + inside the intervals 
X G [b{n — 1), bn\,n = 1,2, ...,N. Here, are the independent limited random quantities with 
a known distribution function P{e), and m < is the negative number sufficiently big to make 
U{x) < at X G [0,A^6]. The length 6 is a specified parameter of the potential U{x). For 
xG[0, A^6], we assume that U{x) = 0. The distribution function P{e) is taken in the following, 
fairly general, form 

P(£) = ^p(^), p{e)>0, Mn = Jp{6)e''de, Mq = 1, Mi = (1) 

The parameter A is the measure of disorder and, at A = 0, the function U{x) represents a 
potential box with a fiat bottom with the depth u and length Nb. At A > 0, we can say that 
the function U{x) corresponds to a potential box with a fiuctuating bottom. Hereafter, we 
imply the tyhermodynamic limit N ^ oo. 

Consider the motion of a particle in such a random potential and formulate the following 
problem. Let the particle, at t = 0, to be located in the point r = Nb, (i.e., at the right side of 
the potential box with the fiuctuating bottom) , and we are seeking for the density of probability 
that the particle will remain in this point at t — )■ oo. From the mathematical viewpoint, it 
means that, at t = 0, the wave function of the particle had the form "^{t = 0,x) = 6{x — r), and 
we have to find D = limt_j.oo(|^(i, The angle brackets here and below indicate averaging 
over realizations of the random potential U{x). This problem is well known in the theory of 
Anderson localization |ll[8] and it can be shown [6] that, if ^„(a;) are the eigenfunctions of the 
Hamiltonian H = —(P/dx'^ +U{x), then D = {J2n \'^n{f)\^)- Nonzero value of D indicates that 
there exist localized functions (i.e., functions essentially nonzero in some finite region, with its 
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size independent of N, at N ^ oo) among eigenfunctions of the Hamiltonian [6l [H [8]. To 
judge about the presence or absence of the locahzed states in the energy interval [U, U + dU], 
in [6] there has been introduced the participation function W{U) defined by the relationship 
W{U)dU = {J2n,Ene[u,u+du] li^nir)]"^) , where En is the eigenenergy of the Hamiltonian H. If 
the eigenfunctions of the Hamiltonian H with energy U are delocalized, then W{U) = 0. 
Otherwise, W{U) is nonzero. For this reason, the participation function will be below referred 
to as spectral dependence of the degree of localization. The main results of this paper are the 
following expressions for the function W and the quantity D: 

These formulas are applicable to the above ID continuous model with a peacewise-constant 
random potential. 

Concluding the introduction, note that the Helmholtz equation describing propagation of 
electromagnetic waves in a layered system, in fact, coincides with the Schrodinger equation 
studied in this paper. This gives the grounds to assert that the results obtained in this paper 
can be used in studies of propagation of electromagnetic waves in ID photonic crystals in the 
presence of disorder. 



2 Continuous model. General properties of the Schredinger 
equation Green's function. 

To solve the above typical problem of the theory of disordered systems, we will apply the 
method of joint statistics of the advanced and retarded edge Green's function (EGF) used in 
[HI E] for analysis of the discrete ID disordered models. A crucial point of the above method is 
the fact that, in the discrete ID model, the EGF of a chain with a single structural unit added 
can be expressed algebraically through the EGF of the initial chain |H|Tl [7] . In this section, we 
will briefly remind the properties of Green's function of the differential Schrodinger equation 
with the Hamiltonian 

H = -d^/dx^ +U{x) (4) 

and, in the next one, we will present a similar relation for the EGF, valid for the case of the 
continuous model with the peacewise-constant potential U{x). 
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Green's function of operator Q is defined by the formula 

G^,^n).j:t^^^^f,/-^^§>^ (5) 

where ipn{(pp) and En{£p) are, respectively, the eigenfunctions and eigenvalues of operator 
Q, corresponding to the discrete (continuous) spectrum. Here, n{p) is the discrete (continuous) 
number of the eigenfunction. The energy argument Vl of the Green's function is, generally, a 
complex number Vt = U — iV (with U and V being real). Using Eq. ([s]), one can show that 
the solution \E'(i(:, x) of the time-dependent Schrodinger equation, with the initial condition 
^(0,x) = ^E'o(a;), can be expressed in terms of the Green's function (|5]), as follows 

*(t,x) = ^lim^ ^ j dUdx'e'^' G^^,{U - iV) ^o{x'), t > (6) 

It can be easily shown that the Green's function (|5| satisfies the differential equation 
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G,,.(n) = *(i-i') (7) 



where only solutions vanishing at x — )■ ±oo should be taken to provide convergence of integrals 
([6]). Taking into accunt the above properties of the Green's function, one can see that the 
quantity D, introduced in Sect.l, is expressed through the product of diagonal elements of the 
advanced and retarded EGF in the following way 

D = lim (|^(t,r)|2) = hm lim / dUM e'^^'-^'^'{Grr{Ui - iVi)Grr{U2 + zV^s)), (8) 

Hear r is the coordinate of the right side of the potential box with a fluctuating bottom 
described in Introduction. Using spectral expansion ([s]), one can show that the quantity D is 
determined only by discrete states of the Hamiltonian (|4]), with the following formula being 
valid [6j 

Z^= lim mt,r)\') = {j:\Mr)\') (9) 

By limiting the region of integration in ([s]) so that Ui,2 E [U,U + dU], we can obtain the 
participation function W{U), introduced in |6] 



W{U)dU=( E \Mr)n (10) 

^ EnG[U,U+dU] ' 

As was already mentioned, nonzero value of W{U) indicates presence of localized states in the 
energy interval [U, U + dU]. In the opposite case, W{U) = 0. 
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3 Case of the peacewise-constant potential. Recurrent 
relations for the EGF. 



Consider the following family of random potentials constant within the intervals of length b: 

( u + En, X E \b(n — l),bn\, x < mb (m<0,?t,- ) , ^ 

U^{x) = \ ^ h U - ^^^^ 

[ U, X > mb 

Here Sn are the independent bounded random quantities with the distribution function P{e) 
([!]), and u is the negative number sufficiently large to meet the condition u+Sn < 0. Assume the 
EGF G'™b „^b(f2) = 7m(^) of the Schrodinger equation (jij) with the potential (2;) to be known. 
Let us pass to the potential Um+i{,x) and consider the EGF G^"^^!) fe(m+i)(^) = 7m+i(^) Eq. 
(4) corresponding to this potential Um+i{x). In this section, we will express the EGF 7m,+i(^) 
through 7m,(^) using the fact that, at x < m6, these potentials are the same. 

Note, first of all, that the discrete spectrum of operator Q with potential (11), we are 



interested in, is positioned on the negative semiaxis. For this reason, in what follows, we 
will consider real part U of the energy argument of the Green's function to be negative VL = 
U — iV,U < 0. The Green's function meets the equations 



[VL + dPldx"^ - W™(x)]G™f,^(r]) = 0, at X < 6m 
[Vl + c^Vc^x2]G™^^(^]) = O, at x > 6m 



(12) 



Let us introduce the functions \I/±(x), so that 

[fi + rfV(^a;2 -W„(x)]^_(x) = 0, with ^_(-cx)) = ^_(m6) = 1 (13) 
[VL + / dx'^]^! +{x) = 0, with ^+(oo) = ^+(6m) = 1 (14) 



Equations and conditions ( 13 ) and ( 14 ) determine the functions (x) in a unique way. It 



follows from Eq. (14) that 

^+(x) = e*^[^-''™l, Q = U-iV, V>0, U <0 (15) 

The Green's function we are interested in can be expressed in terms of the functions (x) 
as follows 

{A^_(x), at X < m6 
(16) 
-B^+(x), at X > m6 

The continuity of the function G^fj^{Q) at x = m6 together with a unit jump of its derivative in 
this point lead to the system of equations for the constants A and B. Solving this system and 
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taking into consideration that \l/-|-(m6) = 1, we obtained for the EGF the following 

relation 



1 



(17) 



Note that the function \E'+(x), entering this relation, is known in the explicit form (15). 

Now, let us pass from the potential Um{x) to the potential Um+i{x) and consider the Green's 
function G^^l^j^^^{^l) . It satisfies the equations similar to (12) 

+ d^dx^ - l{mix)]G2t(ln+i)i^) = 0, at X < 6m 
[n + d'^/dx^ - v]G2Hln+i)i^) = 0, at 6m < X < 6(m + 1), ri = u + Em+i (18) 



[n + d^/dx']G%l^,,{^) = 0, at x > 6(m + 1) 



By analogy with (16), we can write the following expressions for 

A^^ix), at X < hm 



G"!tL = Ce*^^^ + Fe-'^^"", at bm < x < b{m + 1) 



(19) 



G"!tl = B^+ix), at X > 6(m + 1) 



At X = mb, the Green's function G''^h(Tn+i)(^) should be continuous together with its first 
derivative, while, at x = (m + 1)6, the function C™^^^^-) (fi) should be continuous, and its 
derivative should experience a unit jump. This yields four equations for the constants A, G, F, 



and B entering (19). The EGF 7m+i(f2) of interest, corresponding to the potential Um+i{x), 
can be obtained from ( 19 ) 

7„.+i(l^) = G-+^i),,(^-,i)(f^) = B^Mm + 1)) (20) 
Finding the constant B from the above system of equations for A, G, F, and B and using 



Eq. (15) for the function \E'4.(x), we obtain for the EGF 7^+1(1]) the following relation 

/ Q 



ty/n-T]- ^'_(m6) 
y/il - 7] + t (mb) 



+ 



1 



t = tg [bJn - 7]] 



(21) 



Now, using (17) we can express \t''_(m6) through the EGF 7m(fi) 

1 



vl>'__(m6) = iVfi 



(22) 



Here, we took into consideration that ^'^(6m) = ty/Q at U < 0. With the use of Eqs. 



(21) and (22), we can obtain the sought relation between EGF 7m_|_i(fi) and EGF 7m(fi) (the 
corresponding operation will be further referred to as TZ^^) 



7m+l = — — = T^n^r,iln 



(23) 



where 



h 



t 



trj 



t = tg 



VL = U ±lV, ?7 < 0, 1/ = +0, TI = U + Em+l 



bJn-T] 



Below, we will need the operation TZ inverse to (23), which has the following form 

7m+ig - h 



7n 



V7m+1 



m+lj 



(24) 



Thus, the whole method of analysis of the joint statistics of the EGF, developed in [6l [7] for 
discrete models, can be applied to the considered case of a continuous model, corresponding 
to Schrodinger equation Q with the peacewise-constant potential (11). Relevant calculations 
are presented in the following sections. 



4 Calculating spectral dependence of the degree of lo- 
calization 

4.1 Joint statistics of the Green's functions 

Spectral dependence of the degree of localization W{U) and the probability D to find the 
particle at the edge of the random ID system under consideration are given, respectively, by 
Eqs. (10) and ([9]). These formulas are identical to those for similar quantities of the discrete 



modeljniE]. For this reason, for realisation of Eqs. (10) and ([9]), one can use the method based 
on calculation of joint statistics of the advanced and retarded Green's functions developed in 
[HI [7j. Let us remind briefly this method. In the integrals entering Eq. (|8]), we make the 
following substitution u = U2 — Ui,U = Ui and change the notations Grr — ^ 7. Then, for the 
value D (g, we can write the following expression 

D = At lim lim / dudUe'^H-iiU - iVi)-i{U + u + z^)) (25) 

47r^ i-s-oo Vi,2^+0 J 

If the function p{xiyiX2y2) is the joint statistics of the advanced and retarded Green's functions 



entering (25) (here, the arguments Xi and i = 1,2 correspond to real and imagenary parts 



of these functions), then the averaged value of their product can be represented in the form 
{-fill -%Vi)'y{U + u + iV2)) = J dxidyidx2dy2pixiyiX2y2)[xiX2-yiy2 + i{xiy2 + X2yi)] = (26) 

= {X1X2) - {ym) + lixm) + i{yiX2), 
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It was shown in [6] that it suffices to calculate only t{xiy2) and to multiply the result by 4. 
In accordance with [7], calculation of this contribution at Vi^2 +0 can be performed using 
the formula 

{xm) = J dedxidx2P{e)(ru,u+U^iX2)xi{xi)y2ix2), (27) 

with the form of the dependences Xi(xi) and y2{x2) being determined by the fractional-linear 
function 7^^,^(a;) g [] 



h + x ^ an,n + &Q,7?a; 
q + vx~ cn^n + gn,rjX '' 



(28) 



which looks like this 



xi = Re 



2/2 = Im 



, ri = u + e, ill = U - iVi, n2 = U + oj + iV2 (29) 



Function au,u+Lu{xiX2) in Eq. (27) represents the joint statistics of the real Green's functions 



with the energy arguments U and U + u, respectively. In [7J, it has been shown that this 
function satisfies the following equation 



crc/i[/2(xiX2) = J de Pie) 



U2 



'T^Uiri{Xi),TZu2'q{x2) 



dTlu^nixi] 



dxi 



dTZu2-n{^2 



dx-2 



j^ = u + e (30) 



with the operation lZu^r^{x), in this case, being determined by Eq. (24). As shown in [7J, at 



Vi_2 — !■ +0, the relations (29) lead to the following expressions for 2/2(^2) and Xi(xi): 



y2[,X2] 



9u+uj,ri 



5{X2 + 



CU-\ 

gu- 



'Uj,ri 



■u),ri 



Xi(Xi, 



_ au,ri + hu^nX 

yi-s>+0 Cu,-q + gu,riX 



(31) 



(32) 



Substitution of these expressions into (27), allows us to obtain, for the quantity {xiy2) of 
interest the following relationship [7] 

{xiy2) = TT lim / auu+uj{x, a)xdx (33) 

a— !>oo J 

When deriving this relationship, we took into account that the function <7u,u+uj{xiX2) meets 



Eq. (30). It follows from the above formulas that the quantity D, we are interested in, may be 



represented in the form 



D = ^ lim d'^Hxiy2)dudU = - lim e'-'^^a^auu+Jx,a)xdxdudU (34) 

Vr^ Vi,2-!>0,i-s>oo J TT a-s-oo,t-s>oo J 



^We will present this function using notations of [7] 
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As shown in [6j, the participation function W{U) can be obtained from Eq. (34) by omitting 
the integration over U : 



W{U) 



— hm 

7f a— ^oOjt— >oo 



a (7u,u+ujix, a)xdxduj 



(35) 



In this way, the problem is reduced to solving Eq. (30 ). The perturbative approach to equations 



of the type (30 ), proposed in [6j, is the power expansion in A (see Eq. (jTj)), with the first nonzero 
correction being of the order of A^. It was also shown in |j6j that, to calculate the quantities D 
and W{U), only the part of the solution of the equation for the joint statistics (in our case, Eq. 
(30)), singular in co, is needed, with the singularity being of the pole type. Thus, the needed 
singular part (referred to as sing) can be represented in the form 



sing auuj{xiX2) 



A' 



J'u{xiX2) + 0{A') 



(36) 



Now, using Eq. (35), for the function W{U) and quantity we obtain the following 
formulas 



W{U) = -A^ lim / J^u{x, a)xdx + 0(A 



D 



W{U)dU 



(37) 



In the following section, we will present the perturbative approach to Eq. (30) and will derive 



an explicit expression for the function J^u{xiX2) entering Eq. (36). 



4.2 Perturbative approach to Eq. (30) 



To solve the functional equations arising in the perturbation theory described below one has 
to find the eigenfunctions and eigenvalues of the functional operator which acts upon an 

arbitrary function f[x) as follows]^ 



(38) 



where IZ^^j^ is given by (24). We assume the parameter A to be small and represent the sought 
function o'u^u^{xiX2) as a power series in A. 



(JUiU2{^l^2) = QniXi,X2)A'' 
n=0 



(39) 



Let us expand the function auiU2 [T^Ui,ri{xi) ,TZu2,vi^'2)] I '''''^"^^^^^^ \ \ '^'^^2^^^^^ |^ in the right-hand 



side of (30) into a power series in e. Then, Eq. (30) yields 



£g„(xi,a;2)A" 



(40) 



n=0 



^This problem is solved in Appendix, and, in what follows, we will use the results obtained in it. 
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1^ -y 

n,Z=0 



dn 



dxi 



dxo 



e=0 



Remind that the dependence on e in this equation is provided by the quantity t] = u + e. 
By equating the coefficients at the same powers of A in the right- and left-hand sides of Eq. 



(40), we obtain the recurrent relations for the function Qn 



Qo{xiX2) - Qo 



^c/i,n(a;i),7^[/2,n(a;2) 



dTZu^^ujxi] 
dxi 



dTZu2,u{x2) 
dx2 



(41) 



Since the first moment of the function P{e) ([T]) is zero, the quantity Qi vanishes. 



M2 



Q2{XiX2) - Q2 

Qo 



dn. 



dxi 

dnu2,n{x2) 
dx2 



dnu2,u{^2) 



dXn 



(42) 



1] = u + e 



e=0 



and so on. From Eqs. (41 ) and (42 ), we see that they contain the functional operator 'Hui,u, i 



1, 2 (38). Taking into account its properties, described in Appendix, we can immediately write 



the solution of Eq. (41) for Qq: 



Qo{XiX2) — ^Ui,u{Xi)Cir2,u{x2) 



(43) 



To solve Eq. (42), we will present the sought function Q2{xiX2) in the form of expansion over 



eigenfunctions ( 76 ) of operator ( 38 ) : 



Q2{XlX2)= E azS^^'"(Xl)sr^'"(x2) 



(44) 



By substituting this series into the left-hand side of Eq. (42) and by expanding its right-hand 



side using (75), we obtain, for the coefficients Cni, the following formulas: 

Mo 92 ^ 



a 



1 



nl 



l-Xn{Ui,u)Xi{U2,u) 2 de^ 
where the quantities Jn{U e) are defined as 



UUie)Ji{U2e) 



£=0 



(45) 



Jn{Ue) = J 



dx 



dz = J\{Ue), r] = u + e (46) 



Definitions of the functions C and Q entering these expressions are given in Appendix. When 



expanding the right-hand side of (42), we used expression (43) for the function Qo{xiX2). As 



was pointed out above, we are interested only in the part of Q2{xiX2) singular in 00 = U2 — Ui. 



To extract this part, one has to retain, in (44), only the terms with n = —I [B], since only for 



these terms the denominator 1 — Xn{Ui, u)Xi{U2, u) in (45) turns into zero at uj = U2 — Ui = 0. 
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The calculation identical to that performed in [6] leads to the following expression for the 

2 



function J^u{xiX2) entering Eq. (37): 



Jn{Ue] 



£=0 



n 



Now, let us present explicit expressions for the integrals (46): 



Jo(t/,e) = l, uu e) = g^i\nij](f 



r]=u+e 



[Ji(f/,£)]", n>0, f 



Vu -VTT 



u 



lU 



(47) 



(48) 



These expressions are obtained by integrating (46) with the help of residues. Note that, for 



calculations of the derivatives entering Eq. (45), the value e can be considered so small that it 



does not affect positions of the poles of the integrants with respect to the real axis (above or 



below). Using Eq. (77), we can obtain the relationship 



UU) — 



Xn{U,u) 



r — z 



(49) 



which shows that J„(?7, 0) = at n 7^ and that, in the general case, the power expansion of 
Qij]J(^i]]r]<^y) ^ starts from the first power and may be written in the form 



Ji{U,e) = Kue + 0{e' 



(50) 



rj=u+e 



Substitution of this expression into (47) shows that, in sum (47), only the terms with 



J±i{Ue) survive, for which the second derivative of their module squared is nonzero at e = 0. 



Thus, Eq.(47) for the function J^[/(xiX2), can be represented in the form 



-u\Ku? 



U.uf \ U.u/ \ U.uf \ U.uf \ 
Si (Xi)s_'i (X2) - S_'i (Xi)Si' (X2j 



(51) 



Direct algebraic calculations using explicit expressions (75) for function Qu,u{'x) and (23) 
for the operation IZ^j^^i^x) show that 



u 



sm 



u\ 



\Ku? 



sin^fftTfT 



u 



(52) 



2u{U-u) "J' 4(f/-M)2 

Finally, using expressions for the moments and limiting values of the functions a^{x), given in 



[H] (see Appendix), with the aid of Eq. (74), we can obtain the following expressions for the 
first moments and limiting values of the s-functions: 



s^''^{x)xdx = I 



n 



VU -u \t\ 



lim ah^^'^ia) 



\n\ 



u 



71 



U 



t 



(53) 



Then, using Eq. (37), we obtain, for the participation function W{U) and the quantity D, 



expressions ^ and 
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5 Numerical experiment. Localization length. 



The most convincing way to verify theoretical results related to ID solid-state disordered models 
is, nowadays, to compare them with a numerical experiment. Below, we present the results 
of numerical verification of Eq. ^ for the participation function W{U) (spectral dependence 
of the degree of localization), which is considered to be the main result of this paper. In this 



verification, we used definition (10) at dU « \u\. The wave functions entering Eq. (10) 
were obtained by solving numerically the edge problem for Schrodinger equation ^ with the 
random potential (11) using the transfer matrix technique. In the calculations, we assumed 



p{e) = Q{e + 1/2) — Q{e — 1/2) (see Eq. ([T])) and the number of regions of constant potential 
iV ~ 200 - 900. The final function W{U) was obtained by averaging over A^^ ~ 2000 - 4000 
realizations of the random potential. When performing the above calculations, one should 
keep in mind the following: (i) As far as formula ^ obtained in this paper is valid in the 
thermodynamic limit, the number N should be sufficiently large. However, at > 800 — 900, 
in the calculations of the wave functions, the errors arising at multiplications of a great number 
of the transfer matrices rapidly increase; (ii) For a given length bN of a random system, the 
degree of its disorder A, on the one hand, should be large enough for the localization length to 
be smaller than bN, and, on the other, should be small enough not to come out of the range 
of applicability of Eq. (|2]); (iii) In these calculations, one has to check quadratic character of 
the dependence of the computed function W{U) and independence of the results on A^. 

The results of numerical calculations for different values of the parameters b and A are 
presented in Fig.l (the values of all the parameters are given in the figure), the smooth curves 
being calculated using Eq. ^ with no fitting. Figure lb demonstrates better agreement between 
the theory and experiment than Fig. la, because, the above conditions were satisfied much better 
for the case of numerical dependence for Fig. lb. 

In the numerical calculations, it is useful to be able to evaluate the localization length / of 
the wave functions for the random system with a given energy U. For such evaluations, one 
can use the participation function W{U) ^ obtained in this paper. Consider the states of the 
random system with the energies lying within the interval [U,U + dU]. The number of such 
states will be p{U)dU, where p{U) is the density of states. In virtue of spatial uniformity of 
statistical properties of the random potential U{x), we can say that the "centers of gravity" of 
these localized states are distributed more or less uniformly along the x axis. Therefore, the 
number of states in the energy interval [U, U + dU] , whose centers of gravity fall into the spatial 
interval dL of the x-axis, can be estimated as [p{U)/L]dUdL, where L = Nb - is the length 
of the potential box with a fiuctuating bottom. Note that the participation function (10) is 
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mainly contributed by the states whose centers of gravity are separated from the edge of the 
potential box by the distance not exceeding their localization length /. The number dn of such 
states is estimated to be dn = [p{U) / L\ldU . By denoting the mean amplitude of these states 



at a; = as Tpuifi), we can write, for the participation function (10), the following approximate 
expression: 

W{U)dU = |^t/(0)|^dn W{U) = ^|V't/(0)|^ (54) 

With the accuracy acceptable for our purposes, we may assume that the density of states 
p{U) entering this equation does not strongly differ from that po{U) for the potential box with 
the length L and depth u with no disorder. 

p{U) ^ po{U) = £ (55) 

ZTTVt/ — U 



Then, formula (54) yields 



^ViU) - (56, 

The localization length / entering this formula and the amplitude of the wave function ^/'{/(O) 
can be connected by the normalization condition, which will provide a second relationship 
for their calculation. If the quantity ipui'^) were close to a typical amplitude of the wave 
function, the above connection would have a simple form \ipu{{])\^l = 1. However, the arguments 
presented below show that the amplitude of the wave function at x = can be much smaller 
than its typical value, which we denote as ipu. Let us evaluate ipu based on the following 
reasoning. At x < 0, (i.e., outside the potential box), the wave function has the form ipui^) = 



X < I, scattering in the random potential is weak, and the wave 
function, within this interval, approximately corresponds to free motion of the particle with 
the energy U. For this reason, for the wave function near the edge of the random system we 
can write the following expressions 



I ^u{x) = ^u{0) exp[xv^^] X < ^^^^ 
I ^uix) = A sm[x^/U — u + yj] < X < / 

The energy U is assumed here to be sufficiently high, so that the motion of the particle 
has a ballistic, rather than tunnel, character f/ — u — > 0. In the case of small disorder, 
when U — u — 6i^U — u, this requirement does not essentially restrict our consideration. The 
conditions of continuity of the wave function and its first derivative at x = allow one to find 
the constants A and ip: 

(MO)=Asinv ^ (A' = \MO)\'^ 

\tlJu{0)V^ = AVU^ cos ltg<^- ^ 



u 
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The typical values of the wave function module squared in the region of localization /, we 
are interested in, can be estimated as a half of its peak value within [0, /]: 



A 



— maXa.g[o,«] sm^[x\/U - u + ^p] 



(59) 



Using Eq. (58), we can obtain for \'ipu\'^ the following expressions 

l^t/(0)P 



\M0)\' 



u 



u 



when 



IVU -u + io > 



TT 



(60) 



u 



sin^ [ly/U — u + ip] , 



when 



2 u-U 

Now we can apply the normalization condition mentioned above | = 1: 



WU - U + Lp < 



TT 



U-U 



when 



TT 



IVU - u + ip> - 
2 



I^i/(0)p/ u 



u 



U 



sm 



\iVu^ + ^] = 1, 



when 



WU - M + v? < 



TT 



(61) 



(62) 



Equations (61), m2\ and (56) allow us to express the localization length / through the par- 



ticipation function W{U) obtained in this paper. For instance, Eqs. (61) and (56) give the 
following expression for the localization length: 

2|?7-u|3/2 



(63) 



7riy(f/)M2 

For algebraic consistency, we retained here the numerical factor 2/ti. This formula is applicable 
provided that the localization length /, obtained with its aid, meets condition (61): l\/U — u + 



If > ^. Combining Eqs. (62) and (56), we obtain equations for determination of the localization 
length in the case when ly/U — u + y9 < |: 



sm 



u 



arctgi 



lu 



U' 



U 



(64) 



In the topical case when U — u > A, formula (63) appears to be the main one. This is why we 



will not analyze the transcendent equation (64). Figure 2 shows a typical form of wave functions 
of the random system at different energies. Horizontal thick lines show the localization lengths 
obtained using (Q and (|2|. It is seen from Fig. 2 that these formulas may be used to evaluate 
spectral dependence of the localization length for the states of the random system considered 
in this paper. 
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6 Conclusions 



The perturbative approach to the joint statistics of the advanced and retarded Green's func- 
tions, developed previously for the discrete random ID models [6l[7j, is applied to analysis of the 
continuous disordered model described by the Schrodinger equation with a piecewise-constant 
random potential. Using the developed approach, we derived the expression for spectral depen- 
dence of the degree of localization in the sense of the Anderson criterion. Numerical verification 
of the results obtained is presented. In conformity with the commonly accepted opinion, the 
states with negative energies of the considered random system prove to be, generally speak- 
ing, localized, because the participation function (|2|, at these energies, is nonzero. Exceptions 
are the points of delocalization arising at h\/ —u > n (see Eq. ([2|). Unfortunately, we have 
not managed to study behavior of the participation function at large values of the param- 
eter b (when these points appear), because the used algorithm of numerical solution of the 
Schrodinger equation became unstable. In this connection, it makes sense to pay attention to 
similarity between the continuous model described in this paper and the discrete model with a 
complex structural unit [7]. There are strong grounds for believing that the behavior of the par- 
ticipation function W{U) of the continuous model under consideration qualitatively coincides 
with that for the discrete model for which the numerical analysis appears to be feasible. 
In conclusion, emphasize once again that the developed approach and the results obtained can 
be useful for analysis of propagation of the electromagnetic waves in structures of the type of 
ID photonic crystals in the presence of disorder. 



7 Appendix 

Solution of the spectral problem for the operator T-L^i^r] 

To solve this problem, we use the system of eigenfunctions crg.(x) for the operator Hcf{x) = 
f{C — l/x)/x^, obtained in [9] in explicit form. Consider some of these functions cr{x) and 
denote the corresponding eigenvalue of the operator He by A. Then, the following relationship 
should be vahd 

\a{C - 1/x) = Xa{x) (65) 

Let us pass, in this equation, to a new variable y = [x — A]/B, x = A + By, where the 
parameters A and B are supposed to be defined later. If we now introduce a function $(y) = 
C — 1/x = C — + By], then we can easily see that, by passing to the variable y in Eq. 
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( 65 ) , we come to the following relationship 

1 rf$ 



B dy 



a[^{y)]=Xcr{A + By) 



(66) 



Now, we introduce the function s{y) defined as 

s{y) = B a{A + By) and, consequently, a{z) 



1 fz- A 
s 



B \ B 



It follows from (66) that 



1 d$ /^{y)-A 
Bdy ^[ B 



If we define a function TZ by the relation 



Hy)-A 

B 



AB 



1 C-A 1 




r B 1 


+ A H 


/ 


' + 1". 



then Eq. (68) can be rewritten in the form 



dy 



s[n{y)] = Xs{y) 



(67) 



(68) 



(69) 



(70) 



Let us now choose the parameters A, B, and C to make operation (69) coincident with (24). 
This gives rise to a system of equations for these parameters. Solving this system we have 

y/n-ri + tty/Tl tT] „ 2 



A 



B 



C 



t = U 



bJn-ri 



(71) 



f(n-r]){l + t^) ^(fi-r/)(l + t2) ^ Vl + t^' 

Thus, the function g(y)([67)) constructed with the aid of the eigenfunction a{x) of the operator 



He, for the parameters A, B, and C determined by Eq. (71 ) is the eigenfunction of the operator 



^Q,r; (38), with the appropriate eigenvalue A being coincident with that of the operator He- 
Below, we present a compact expression for the s-functions. 

As shown in |:9], an arbitrary function f{x) may be expanded in series in terms of the 
functions erg- (a;). Remind the explicit form of the functions cr^(a;) and expressions for the 
eigenvalues A„ of the operator Hq at C < 2: 



Lc{x)G''{x), Lc{x) = 



2tii 



X — R X — R* 



, G{x)^ 



R* 



X 



R-x 



(72) 



A. 



R 



R* = 



C-tV4:- C2 



\C\ < 2 



as well as the rules of expansion of an arbitrary function f{x) in series in terms of the above 
functions: 



+00 

f{x)= J2 ^naS(x), 



G^(x) 



-dx 



(73) 
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Using these relationships, we can obtain similar rules for the eigenfunctions (67) of the 



operator(38) (called by s-functions) 



^n.^s^-n?/) = A„(fi, V)s^''iy), s^'^iy) = B al{A + By), A„(fi, r,) = exp[2tnb^Q - v] 

(74) 



Here, the parameters A, B, and C are defined by formulas (71). The superscript of the s- 



functions indicates their dependence on the energy argument Q. 
The arbitrary function f{x) can be represented as the series 

' f{y) 



{A + By) 



dy 



fiy) 



-dy, 



with 



r ( ^- ^* -y u - V^-V^-V , 
yn,r,[y) == , where r = r 



Gn,r,iy) 



(75) 



IT] 



r — y irj 
Using the quantities introduced in this way, we may write the following compact expressions 
for the s^''^{y) s-functions: 



s^'^(y) = CnM0n,r,iy), ^nM = 



2m 



y — r y — r 



(76) 



Here, the subscripts of the Lorentzian C indicate dependence of this function on the energy 
parameters Q and 1] = u + e. One can easily make sure that 



Sn,,iT^n,r,ix)) = X^{n,rj)g-^^{x) 



(77) 
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Figure 1: Spectral dependence of the degree of localization of the states for the ID disordered 
system with a piecewise-constant random potential. The noisy curves are obtained by computer 
simulation and the smooth ones are computed using Eq. 
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Figure 2: Wave function of the ID disordered system with a piecewise-constant random 
potential for different values of the energy U. Thick horizontal lines show the localization 



length calculated using Eq. (63). 
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